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Abstract. To demonstrate the role played by the von Neumann entropy spectra 
in quantum phase transitions we investigate the one-dimensional anisotropic 
S\J{2)0XXZ spin-orbital model with negative exchange parameter. In the case 
of classical Ising orbital interactions we discover an unexpected novel phase 
with Majumdar-Ghosh-like spin-singlet dimer correlations triggered by spin-orbital 
entanglement and having fc = 7r/2 orbital correlations, while all the other phases are 
disentangled. For anisotropic XXZ orbital interactions both spin-orbital entanglement 
and spin-dimer correlations extend to the antiferro-spin/alternating-orbital phase. 
This quantum phase provides a unique example of two coupled order parameters which 
change the character of the phase transition from first-order to continuous. Hereby we 
have established the von Neumann entropy spectral function as a valuable tool to 
identify the change of ground state degeneracies and of the spin-orbital entanglement 
of elementary excitations in quantum phase transitions. 
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1. Spin-orbital physics and von Neumann entropy spectra 

In the Mott-insulating limit of a transition metal oxide the low-energy physics can be 
described by Kugel-Khomskii-type models [1], where both spin and orbital degrees of 
freedom undergo joint quantum fluctuations and novel types of spin-orbital order [2] 
or disorder [3] may emerge. Following the microscopic derivation from the multiorbital 
Hubbard model, the generic structure of spin-orbital superexchange takes the form of a 
generalized Heisenberg model nig, 

Y. , (i) 

(*t>ll7 

as indeed found not only for the simplest systems with S' = 1/2 spins: KCuFs [1], the 
RTiOs perovskites [6], LiNi02 and NaNi02 [Z], Sr2Cu03 [8], or alkali RO2 hyperoxides 
[9], but also for larger spins as e.g. for S' = 2 in LaMnOa [10]. In such models 
the parameters that determine the spin-S' Heisenberg interactions stem from orbital 
operators and — they depend on the bond direction and are controlled by 
the orbital degree of freedom which is described by pseudospin operators {Tj}. That 
is, these parameters are not necessarily hxed by rigid orbital order Bini, but quantum 
fluctuations of orbital occupation [I21[I3] may strongly influence the form of the orbital 
operators, particularly in states with spin-orbital entanglement (SOE) [TH[T5]. As a 
consequence, amplitudes and even the signs of the effective exchange can fluctuate in 
time. Such entangled spin-orbital degrees of freedom can form new states of matter, 
as for instance the orbital-Peierls state observed at hnite temperature in YVO3 [T61IT7] . 
Another example are the collective spin and orbital excitations in a one-dimensional 
(ID) spin-orbital chain under a crystal held which can be universally described by 
fractionalized fermions na. It is challenging to ask which measure of SOE would be 
the most appropriate one to investigate quantum phase transitions in such systems. 

The subject is rather general and it has become clear that entanglement and other 
concepts from quantum information provide a useful perspective for the understanding 
of electronic matter [T91123] . Other examples of entangled systems are: topologically 
nontrivial states m, relativistic Mott insulators with 5d ions ESI, ultracold alkaline- 
earth atoms [26], and skyrmion lattices in the chiral metal MnSi [27] . 

One well-known characterization of a quantum system is the entanglement 
entropy (EE) determined by bipartitioning a system into A and B subsystems. This 
subdivision can refer for example to space [I9l[28], momentum [28l|29], or different 
degrees of freedom such as spin and orbital PI. A standard measure is the von 
Neumann entropy (vNE), = —Ai A{p\^og 2 p\}, for the ground state |To) which 
is obtained by integrating the density matrix, p\ = Trs|To)(To|, over subsystem B. 
Another important measure is the entanglement spectrum (ES) introduced by Li and 
Haldane [3T], which has been explored for gapped ID spin systems [32], quantum 
Heisenberg ladders [33], topological insulators [3l], bilayers and spin-orbital systems [28] . 
The ES is a property of the ground state and basically represents the eigenvalues Pi of 
the reduced density matrix p\ obtained by bipartitioning of the system. Interestingly 
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a correspondence of the ES and the tower of excitations relevant for SU(2) symmetry 
breaking has been pointed out recently [35]. It was also noted that the ES can exhibit 
singular changes, although the system remains in the same phase [36] . This suggests 
that the ES has less universal character than initially assumed [37] . 

In this paper we explore a different entanglement measure, namely the vNE 
spectrum which monitors the vNE of ground and excited states of the system, for 
instance of a spin-orbital system as dehned in equation (1). In this case we consider 
the entanglement obtained from the bipartitioning into spin and orbital degrees of 
freedom in the entire system jSO]. Here the vNE is obtained from the density matrix, 

= Tro|T„) (Tji|, by taking the trace over the orbital degrees of freedom (Tro) for 
each eigenstate We show below that the vNE spectrum, 

5.n(w) = -Yi log2 {‘^ - g'n} . (2) 

n 

reflects the changes of SOE entropy for the different states at phase transitions. The 
excitation energies, Un = En — Eq, of eigenstates |\k„) are measured with respect to 
the ground state energy Eq. It has already been shown that the vNE spectra uncover 
a surprisingly large variation of entanglement within elementary excitations ffl. Also 
certain spectral functions have been proposed, that can be determined by resonant 
inelastic x-ray scattering [38|, and provide a measure of the vNE spectral function. 
Here we generalize this function to arbitrary excitations |d^n), i-e., beyond elementary 
excitations which refer to a particular ground state. We demonstrate that focusing on 
general excited states opens up a new perspective that sheds light on quantum phase 
transitions and the entanglement in spin-orbital systems. 

The paper is organized as follows: In section [2] we introduce the ID spin-orbital 
model with ferromagnetic exchange, and in section [3] we present its phase diagrams 
for the Ising limit of orbital interactions and for the anisotropic S\J{2)^XXZ model 
with enhanced Ising component. SOE is analyzed in section 0] using both the spin- 
orbital correlation function and the entanglement entropy and we show that these two 
measures are equivalent. In section [5] we present the entanglement spectra and discuss 
their relation to the quantum phase transitions. The main conclusions and summary are 
given in section | 6 l The distance dependence of spin correlations in the antiferromagnetic 
phase is explored in the Appendix. 

2. Ferromagnetic S\J{2)<ZXXZ spin-orbital model 

The motivation for our theoretical discussion of spin-orbital physics comes from t 2 g 
electron systems in which orbital quantum fluctuations are enhanced by an intrinsic 
reduction of the dimensionality of the electronic structure na. Examples of strongly 
entangled quasi-lD t 2 g spin-orbital systems due to dimensional reduction are well known 
and we mention here just LaTiOa | 6 ], LaVOs and YVO 3 na, where the latter two 
involve {yz, zx} orbitals along the c cubic axis; as well as p^ and Py orbital systems in 
ID fermionic optical lattices [391 - ITT] . This motivates us to consider the ID spin-orbital 
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model for S = 1/2 spins and T = 1/2 orbitals with anisotropic XXZ interaction, i.e., 
with reduced quantum fluctuation part in orbital interactions. The |+) and |—) orbital 
states are a local basis at each site and play a role of yz and zx states in t 2 g systems, 

L 

'H{x,y,A) = - j5^(4-4+i + a;)([f,-f,+i]A + 2/), (3) 

i=i 

[f, ■ h+ilA = A + TfTJ^,) + T-Tp„ (4) 

where J > 0 and we use periodic boundary conditions for a ring of L sites, i.e., L + 1 = 1. 
The parameters of this model are {x,y} and A. At x = y = 1/A and A = 1 the model 
has SU(4) symmetry. Hund’s exchange coupling does not only modify x and y but also 
leads to the XXZ anisotropy (A < 1), a typical feature of the orbital sector in real 
materials [5l[T2]. The antiferromagnetic model (J = —1) is Bethe-Ansatz integrable at 
the SU(4) symmetric point jl2] and its phase diagram is well established by numerical 
studies [IHlIll] . It includes two phases with dimer correlations [15] which arise near the 
SU(4) point. Some of its ground states could be even determined exactly at selected 
{x,y,A) points [461 - I50] . 

Here we are interested in the complementary and less explored model with negative 
(ferromagnetic) coupling (J = 1 ), possibly realized in multi-well optical lattices [5T], 
which has been studied so far only for SU(2) orbital interaction (A = 1) [30]. This 
model is physically distinct from the antiferromagnetic (J = —1) model, except for the 
Ising limit (A = 0) where the two models can be mapped onto each other, but none 
was investigated so far. The phase diagrams for J = 1, see figure H] determined using 
the fidelity susceptibility [52] display a simple rule that the vNE (|2]) vanishes for exact 
ground states of rings of length L which can be written as products of spin (It/s)) and 
orbital (iV'o)) part, |To) = I'ljjs) 0 

3. Phase transitions in the spin-orbital model 

To understand the role played by the SOE in the ID spin-orbital model ([3]) and (|4]) we 
consider the phase diagrams for A = 0 and A = 0.5, see figure [H In the case A = 0 
all trivial combinations of ferro (F) and antiferro (A) spin-orbital phases labeled I-IV 
have = 0, i.e., spins and orbitals disentangle in all these ground states: FS/FO, 
AS/FO, AS/AO, FS/AO. If both subsystems exhibit quantum fluctuations, the ground 
state I To) can no longer be written in the product form. This occurs for the AS/AO 
phase III at A > 0. 

Most remarkable is the strongly entangled phase V at A = 0 and y < 0, see figure 
[l](a). This phase occurs near x ~ — (5'j• 5'j+i)AF = ln2 — 1/4, i.e., when the uniform 
antiferromagnetic spin correlations in phase III are compensated by the parameter x, so 
that the energy associated with Hamiltonian ([3]) de facto disappears. This triggers state 
V with strong SOE (see below) as the only option for the system to gain substantial 
energy in this parameter range by nonuniform spin-orbital correlations. The analysis of 
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X X 


Figure 1. Phase diagrams of the spin-orbital model [equation Q] obtained by two 
methods, fidelity susceptibility or an exact diagonalization of an L = 8 site model, for; 
(a) A = 0, and (b) A = 0.5. The spin-orbital correlations in phases I-IV correspond to 
FS/FO, AS/FO, AS/AO, FS/AO order (see text). At A = 0 only the ground state of 
a novel phase V has finite EE, > 0 (shaded), whereas at A > 0 the EE in phases 
III and VI is also finite. 


phase V in 


terms of the longitudinal equal-time spin (orbital) structure factors 

1 L 

= 75 E (o;:.oa, 


m,n=l 


(5) 


where O = S or T, reveals in £gure[2](a) at A = 0 and y = —1/4 for the spin structure 
factor S^^{k) oc (1 — cos/c). This is a manifestation of nearest neighbour correlations, 
while further neighbour spin correlations vanish and moreover we hnd a quadrupling in 
the orbital sector, see £gure[2](b). Thus the spin correlations indicate either a short-range 
spin liquid or a translational invariant dimer state. 

The hidden spin-dimer order [53] can be detected by the four-spin correlator (we 
use periodic boundary conditions). 


flW = 7E 


i^j' ^j+l){^j+r' Sj+l+r 


SrSi+,y 


( 6 ) 


j=l 

At A = 0 we hnd \D{r)\ with long-range dimer correlations in phase V, but not in III. 
Phase III is a state with alternating [k = tt) spin (orbital) correlations in the range 
X < 0.17 shown in figures [2](a,b). Interestingly for A > 0 the dimer spin correlations 
\D{r)\ are not only present in phase V but also appear in phase III. Moreover a phase 
VI emerges, complementary to phase V, with interchanged role of spins and orbitals, see 
£gure[I](b). The order parameters for phase VI follow from the form of structure factors 
which develop similar but complementary momentum dependence to that for phase V 
seen in £gure[2](a), i.e., maxima at 7r/2 for S^^{k) and at vr for T^^{k). We remark that 
phases V and VI are unexpected and they were overlooked before for the SU(2)(8)SU(2) 
model at A = 1 [30|. From the size dependence of |71(r)| in hgures [2](c,d) we conclude 
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Figure 2. Top— Spin and orbital [T^^{k)] structure factors ([5]) for the ID 

spin-orbital model ([3|) of L = 8 sites at A = 0 and y = —1/4: (a) S^^(k) and (b) 
T^^{k). Bottom— Spin dimer correlations D{r) equation dH) found at A = 0.5 for 
decreasing 1/r for clusters of (c) L = 12 and (d) i = 16 sites. 


that the dimer correlations are long-ranged at A = 0.5 in phase V, but also in III, as 
seen from the data for x G [0.0, 0.4), where they coexist with the AS correlations. 

These results suggest that the ground state V in figure [I][a) is formed by spin-singlet 
product states 

|<l.f) = [l,2][3,4][5,6]...[L-l,L], 

|•1>?) = |2,3]|4,5]|6,7]...|L,1|, (7) 

where [/, /-b 1] = (| ti) — | it))/\/2 denotes a spin singlet. They are not coupled to orbital 
singlets on alternating bonds as it happens for the AFantiferromagnetic SU(2)(8)SU(2) 
spin-orbital chain in a different parameter regime [16] , but to Ising configurations in the 
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orbital sector. The four-fold {k = tt/2) periodicity of orbital correlations is consistent 
with four orbital states: 


^1) = 1+ + -- + + •• 



+ 

1 

1 

+ 

+ 

1 

•• + -), 


1 
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1 
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1 
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+ 

1 

1 

+ 


(8) 


The decoupling of singlets is complete for ?/ = —1/4 and A = 0, where (++) and (-) 

bonds yield vanishing coupling in equation (E]), and the phase boundaries of region V 
are Xjjjy = 3/4 -|- 2{Sj-Sj^i)x¥ — 0.136 and Xyjj = 3/4 in the thermodynamic limit; 
moreover we hnd perfect long-range order of spin singlets, i.e., D{r) = (3/8)^(—1)^. 

The dimerized spin-singlet state at A = 0 has the same spin structure as the 
Majumdar-Ghosh (MG) state [5l], however its origin is different. While the MG state 
in a J 1 -J 2 Heisenberg chain results from frustration of antiferromagnetic exchange (at 
J 2 = Ji/2), here the spin singlets are induced by the SOE. At A = 0 the only phase with 
hnite SOE 5°^ = 1 is phase V, see hgure El In contrast, for A > 0 one hnds hnite EE 
also in phase III, when the original product ground state changes into a more complex 
superposition of states and joint spin-orbital fluctuations na appear. These correlations 
control the SOE and give equivalent information to see section |H Furthermore, EE 
increases with x towards phase V where it is further amplihed and exceeds = 1. The 
related softening of orbital order will be discussed below. Interestingly we hnd a one- 
to-one correspondence of hnite EE and long-range order in the spin dimer correlations 



X 


Figure 3. Spin-orbital entanglement entropy in the ground state of the spin- 
orbital model © for the three phases III, V and II as a function of x for various 
A. Solid line for A = 0 stands for the k = 0 ground state in the limit of A —0. 
Parameters: y = —0.5 and L = 8. 
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\D{r)\. 

The superstructure of phase V emerges from the interplay of spin and orbitals, 
where orbitals modulate the interaction of spins in equation ([3]), and vice versa. It is 
important to distinguish this from the Peierls effect, where the coupling to the lattice 
is an essential mechanism. The orbital Peierls effect observed in vanadates [lElin] or 
the orbital-selective Peierls transition studied recently [55] fall into the former category, 
yet, as they involve orbital singlets — they are distinct from the case discussed here. 

4. Spin-orbital entanglement 

The description of spin-orbital entanglement in terms of the vNE entropy, as discussed in 
section [3l is a very convenient measure of entanglement. But it is also a highly abstract 
measure. To capture its meaning, one has to refer to mathematical intuition, namely to 
the fact that any product state, |\1/) = I'lps) ® I'lpo), has zero vNE. That is, an entangled 
state is a state that cannot be written as a single product. A more physical measure 
are obviously spin-orbital correlation functions relative to their mean-held value na. 
Such correlation functions vanish for product states where mean-held factorization of 
the relevant product is exact, i.e., spins and orbitals are disentangled. 

To detect spin-orbital entanglement in the ground state we evaluate here the joint 
spin-orbital bond correlation function Ci for the S\J{2)<^XXZ model ([3|), dehned as 
follows for a nearest neighbour bond {i,i -|- 1) in the ring of length L [H], 

Ci = - {((^* • - ( a , • ( 9 ) 

2 = 1 

The conventional intersite spin- and orbital correlation functions are: 

1 ^ 

= ( 10 ) 

2 = 1 

1 ^ 

Tr^-Y,{T^■f,+r). (11) 

2 = 1 

The above general expressions imply averaging over the exact (translational invariant) 
ground state found from Lanczos diagonalization of a ring. While Sr and correlations 
indicate the tendency towards particular spin and orbital order, Ci quantihes the spin- 
orbital entanglement — if Ci 7 ^ 0 spin and orbital degrees of freedom are entangled 
and the mean-held decoupling in equation ([3|) cannot be applied as it generates 
nncontrollable errors. 

Figures |4j(a) and|4j(b) show the nearest neighbour correlation functions Si, Ti and 
Cl sX y = —0.5, for A = 0 and A = 0.5, respectively, as functions of x. The nearest 
neighbour spin correlation function Si is antiferromagnetic (negative) in all phases III, 
V and II shown in hgure H] while {negative) Ti indicates AO correlations in phase III 
and ferro-orbital {positive) in phase II. Finite A = 0.5 triggers orbital huctnations which 
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Figure 4. Nearest neighbour spin CHD, orbital Ti (ED, and joint spin-orbital 
Cl (0 correlations as obtained for a spin-orbital ring (ED with L = 12 sites and 
y = —0.5, as functions of x for: (a) A = 0, and (b) A = 0.5. The III-V phase 
boundary (dotted vertical line) in (b) has been determined by the maximum of the 
fidelity susceptibility [52] . 


lower Ti below the classical value of 0.25 found at A = 0. In the intermediate spin dimer 
phase Ti is negative for all A > 0, while it is zero for A = 0. 

It is surprising that Ci is positive in phase V at A = 0 in spite of the classical Ising 
orbital interactions, see £gure|D(a). It is also positive in phases III and V at A = 0.5 [see 
figure HDb)]. Note that positive Ci is found in the present spin-orbital chain with J > 0, 
while Cl is negative when J < 0 [12]. In phase II Ci vanishes in the entire parameter 
range as then the ground state can be written as a product. The same is true for phase 
III at A = 0. We emphasize that the dependence of Ci on x is completely analogous to 
that of the von Neumann entropy in hgure [3l which also displays a broad maximum in 
the vicinity of the III-V phase transition at A = 0.5, and a step-like structure in phase V 
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at A = 0. Thus we conclude here that the vNE yields a faithful measure of SOE in the 
ground state that is qualitatively equivalent to the more direct entanglement measure 
via the spin-orbital correlation function Ci na. 

5. Entanglement spectra and quantum phase transitions 

Figure [3] stimulates the question about the origin and the understanding of the sudden 
or gradual EE changes at phase transitions. This can be resolved by exploring the 
vNE spectral function dehned in equation ([2]) and shown in hgures [5]^a) and [5]^b) 
for A = 0 and 0.5, where colors encode the vNE of states. The excitation energies 
^n{,x) = En{x) — Eq{x) are plotted here as function of the parameter x. Only the lowest 
excitations are shown that are relevant for the phase transitions and the low-temperature 
physics. They include: (i) the elementary excitations of the respective ground state, 
and {ii) the many-body excited states that are relevant for the phase transition(s) and 
may become ground states or elementary excitations in neighbouring phases when the 
parameter x is varied. 

The AS/FO ground state of phase II in hgure [5|^a) obtained for a ring of L = 8 
sites is an AS singlet {S = 0) with a maximal orbital quantum number, T = L /2 = 4, 
and a twofold {k = 0, tt) degeneracy at A = 0. The spin excitation spectrum appears 
as horizontal (red) lines and consists of gapless triplet S' = 1 excitations. The low- 
lying excitations of the Bethe-Ansatz-solvable antiferromagnetic Heisenberg chain form 
a two-spinon (s-s) continuum, whose lower bound is given by e{k) = 7r|sin/c|/2 in the 
thermodynamic limit [56]. For the L = 8 ring the spectrum is discrete with a Ak = vr/d 
spacing, and it is known that the energy of triplet excitations £‘^(vr) will scale to zero 
as 1/L [57ll5^ . Red lines in phase II with finite slope are orbital excitations. The 
x-dependence is due to the spin part of "H (I3|) which determines both the orbital energy 
scale and the dispersion, J-p = (x-\- {Sj ■ S'^+i) af) (1 — A cos k). This energy changes with 
X and at hnite A also with momentum k, see hgure [5](b). While the orbitons are gapped, 
the low-lying excitations are either magnons or x-dependent spin-orbital excitations. It 
is remarkable that the latter are entangled in general, although the ground state II is 
disentangled. 

With decreasing x a hrst-order phase transition from II to V occurs by level crossing 
of disentangled (red) and entangled (green) ground states. The spin-singlet (S' = 0) 
ground state of phase V has degeneracy 4 at A = 0, and its components are labeled 
by the momenta k = 0, ± 7 r/ 2 , 7 r. This is rehected by hnite 9 ?^ order parameter in 
hgure [ 6 ](a). Note that at A > 0 this four-fold ground state degeneracy is lifted. In the 
spin-dimer phase a gap opens in the spectrum of elementary spin excitations [ 6 n[[ 6 T]. 
The one-magnon triplet gap A 5 (( 5 ) oc 5^/^ depends on y via the dimerization parameter 
5 = 1/|4|/|. In phase V even the pure magnetic excitations are entangled [see horizontal 
green lines in hgure[5](a)]. The lowest excitations in the vicinity of the phase transitions 
have orbital character. From hnite EE in hgure |5][a) one recognizes that these states 
are inseparable spin-orbital excitations. 
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Figure 5. vNE-spectrum of lowest energies En{x) (relative to the ground state energy 
Eq{x)) versus x with colors representing the size of the vNE of individual states. Data 
for the three phases III, V and II is shown for y = —0.5, L = 8 and: (a) A = 0, and 
(b) A = 0.5. Here e®(fc) denotes spin (orbital) excitation, e^{k) corresponds 

to an elementary excitation having the same S and T as the ground state, and e^^(k) 
stands for the spin excitation under simultaneous flipping of orbitals. 


The phase transition from the dimer phase V to the AS phase III (S' = 0) appears 
singular in the sense that it is hrst order at A = 0 and continuous otherwise [hgures 
|5](a,b)]. To locate the center of the continuous phase transition between phases III and V 
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Figure 6. Ground state energy relative to phase V, Eo{x) — Eq (x) (dots), orbital order 
parameters (dashed), -0^ = [T^^( 7 r)] 2 , = [T^^( 7 r/ 2 )] 2 , and bond spin correlations 

I'S'il = |(.Si-.S 2 )| (solid line), for phases III, V and II (from left to right), for; (a) A = 0.0 
and (b) A = 0.5. Parameters: y = —0.5 and L = 12. 


at A > 0, we have selected the peak of the hrst derivative of the entanglement entropy, 
see figure |3l Yet also the peaks in the derivatives of the fidelity susceptibility, the orbital 
correlation function Ti [see figure ll][b)] and the orbital order parameters "0^ and in 
figure |6](b) may be used. Finally, we note that the scaling of entanglement with system 
size has quite different behaviour in phases III and V, indicating that a phase transition 
separates them. 

Furthermore, the peculiar feature of the AS/AO phase III manifests itself in a 
twofold degeneracy and zero SOE at A = 0 in contrast to the nondegenerate ground 
state and finite SOE at finite A. The entanglement has two sources, namely: {i) the 
interplay of quantum fluctuations in the spin and orbital sectors and {ii) the dimerization 
order which coexists with antiferromagnetic spin correlations in phase III at finite A. 
The latter is the origin of the nondegenerate ground state as it yields a coupling to 
the £®°(7r) excitation (nearly horizontal in x), and the emergence of the spin-dimer 
correlations D{r) leads to a faster decay of the spin correlations in phase III than in the 
ID antiferromagnetic Heisenberg chain, see the Appendix. The orbital order parameters 
and compete in phases III and V, see figure |6](b), near the phase boundary in 
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figure [U^b). This also explains why the transition from phase V to III is smooth at hnite 
A in terms of both the vNE (£gure[3]) and the nearest neighbour spin correlations |S'i|. 

6. Conclusions and summary 

Summarizing, we have studied the quantum phases and the spin-orbital entanglement of 
the ID ferromagnetic SU{2)^XXZ model by means of the Lanczos method. We have 
discovered a previously unknown translational invariant phase V with long-range spin 
singlet order and four-fold periodicity in the orbital sector. Its mechanism is distinct 
from the dimer phases found in the ID antiferromagnetic spin-orbital model near the 
SU(4) symmetric point jl5]. Both III-V and II-V phase transitions arise from the 
spin-orbital entanglement in the case of Ising orbital interactions. When the orbital 
interactions change from Ising to anisotropic XXZ-tjpe, the entanglement develops 
in phase III, where antiferromagnetic spin correlations and long-range spin dimer 
order coexist, changing the quantum phase transition from hrst-order to continuous. 
Furthermore in the regime of hnite orbital huctuations (A > 0) another phase VI 
emerges, which is complementary in many aspects to phase V, but with the important 
difference that phase VI disappears in the limit A = 0. 

We have shown that the von Neumann entropy spectral function iSvn(i^) (E]) is 
a valuable tool that captures the spin-orbital entanglement SOE of excitations and 
explains the origin of the entanglement entropy change at a phase transition. From the 
perspective of spin-orbital entanglement we encounter [i) hrst-order transitions between 
disentangled (II) and entangled (V) phases, {ii) a continuous transition involving two 
competing order parameters between two entangled phases. III and V, and {in) trivial 
hrst-order transitions between two disentangled phases. Case {ii) goes beyond the 
commonly accepted paradigm of a single order parameter to characterize a quantum 
phase. 

Moreover, we have presented two simple measures of entanglement in the ground 
state and shown that they are basically equivalent — the direct measure via the (quartic) 
spin-orbital bond correlation function Ci ([9]) and the von Neumann entropy The 
latter is defined by separating globally spin from orbital degrees of freedom in the ground 
state. 
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Figure 7. Modulus of spin correlations Sr equation OT versus the inverse distance 
1/r as obtained for a spin-orbital ring ([3|) with L = 16 sites, for: (a) A = 0, and (b) 
A = 0.5. Parameter: y = —0.5. 


Appendix: Distance dependence of the antiferromagnetic spin correlations 

Here we explore in more detail the competition of the antiferromagnetic (AF) spin 
correlations of the spin-orbital chain in the AS/AO phase III and the Majumdar-Ghosh 
like spin-singlet dimer correlations that coexist at hnite A, as we found in our work. 
For A = 0 the spin correlations in phase III are those of an AF Heisenberg spin chain, 

(s. ■ S.+,) ~ (-ir (12) 

which reveal the typical 1 /r-power law decay combined with logarithmic corrections that 
were hrst predicted by conformal held theory [6^163] as well as by renormalization group 
methods |6l], and subsequently conhrmed j65] by numerical density matrix method [66] . 

In hgure O^a) we present our numerical data for the spin-correlation function Sr 
equation oa (i.e., for translational invariant ground states) for several values of x, 
and for A = 0 and y = —0.5. In the A = 0 case there are only two distinct types of 
behaviour of namely exponential decay in phase V and the power law decay of the 
ID quantum Neel spin liquid state, which are the same in phases H and HI. 

Figure [71)b) displays S'r at A = 0.5 for different x-values. Here again the 
unperturbed AF correlations of the ID Neel spin-liquid state appear in phase H 
{x > 0.7). It is evident that in phase HI the AF spin correlations are strongly reduced, 
due to the competition with the coexisting long-range ordered spin-singlet correlations. 
The spin singlet order increases with x in phase HI, and as a consequence we observe 
here that the decay of Sr becomes stronger as x approaches the HI/V transition. 

In hgure [8] we present a logarithmic plot which highlights the different decays of Sr 
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Figure 8. Logarithm of modulus of spin correlations Sr equation OT for increasing 
distance r as obtained for the spin-orbital model equation ([S]) on a ring of L = 16 sites 
for A = 0.5, y = —0.5, and three values of x. Exponential decay of Sr with increasing 
r is obtained for phase V (x = 0.6). 


for A = 0.5 in the three different phases: III, V, and II. We have selected the values for 
X = 0.3, 0.6 and 0.7, respectively, for greater transparency. The log-plot shows clearly 
the exponential decay of Sr in phase V. It also shows that the L = 16 system reveals 
strong hnite size effects in phase II where Sr has power law decay. Nevertheless it is 
clear already from the L = 16 data that the AF spin correlations in phase III (here 
shown for x = 0.3) are strongly suppressed and approach the exponential decay of Sr in 
phase V (x = 0.5 and 0.6) when x approaches the III-V phase boundary from the left. 

Summarizing, we hnd that in phase III the AF spin correlations of the ID Neel 
spin liquid state decay much more rapidly as the competing spin-singlet order emerges. 
This effect is particularly strong near the boundary of phase III to the spin-singlet 
dimer phase V. Whether in the thermodynamic limit the correlations Sr also decay 
exponentially in phase III as in V cannot be decided here, and this question is beyond 
the scope of the present work. 
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